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ABSTRACT 

The growth of density perturbations in an expanding universe in the non-linear regime 
is investigated. The underlying equations of motion are cast in a suggestive form, and 
motivate a conjecture that the scaled pair velocity, h{a,x) = —[v/{ax)] depends on the 
expansion factor a and comoving coordinate x only through the density contrast a{a,x). 
This leads to the result that the true, non-linear, density contrast < {5p/ p)^. >^/^= a{a, x) 
is a universal function of the density contrast aL(a, /), computed in the linear theory and 
evaluated at a scale / where / = x(l -|- a^)^/"^. This universality is supported by existing 
numerical simulations with scale-invariant initial conditions having different power laws. 
We discuss a physically motivated ansatz h{a,x) = /i[a^(a, x)] and use it to compute 
the non-linear density contrast at any given scale analytically. This provides a promising 
method for analysing the non-linear evolution of density perturbations in the universe and 
for interpreting numerical simulations. 
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1. Introduction 

It is generally believed that large scale structures in the universe formed through the 
growth of small inhomogeneities via gravitational instability. One convenient measure of 
the inhomogeneities present in the universe at any time t is provided by the mean 
square fluctuation in the mass, a'^{t,x) —< {6M/M)1 >, contained within a sphere of 
comoving radius x. When cr^ <C 1, the evolution of inhomogeneities can be studied using 
linear perturbation theory in an expanding background. But when cr^ grows to a value 
of order unity, linear theory breaks down. This phase of the evolution is usually handled 
by numerical simulations and analytic methods seem to have very limited validity. Our 
understanding of structure formation will increase significantly if the results of numerical 
simulations can be understood in the framework of simple physical concepts and analytic 
approximations. This paper discusses one possible approach to understand the evolution 
of density fluctuation in the nonlinear domain. 

The idea behind this approach is as follows: It is known (Peebles 1980) that the 
differential equation governing the correlation function ^{t, x) is essentially 'driven' by the 
average pair velocity fpair(^, x). We present arguments which motivate the ansatz that the 
ratio [vpairit, x)/ — (ax)] depends on t and x only through the density contrast a'^{t,x); 
i.e., we can set Vpair{t,x) — [—axh{a'^)] where h{a'^) is some specific function which is 
independent of the initial conditions chosen. With this physically motivated assumption, 
it is possible to derive a differential equation for and integrate it in closed form. The 
solution shows that, a'^{t,x) can be related to the density contrast a\{t,l) evaluated by 
linear theory at the scale I = a;(l + cr^)^/^ in a universal fashion. Such universality, if 
established, would be a key feature of gravitational dynamics in an expanding Universe. 
Using this result one can compute the density contrast at highly non-linear epochs knowing 
the linear density contrast evolved to the same epoch. 

The paper is organized as follows: In the next section we discuss Newtonian mechan- 
ics in an expanding universe and stress certain properties of the resulting equations. In 
particular, as is well known, they admit the possibility of similarity solutions in which the 
statistical properties of density fiuctuations at different times are identical apart from a 
change of spatial scale. In section 3 we discuss the similarity solution in some more detail 
and motivate possible generalizations of this principle to estimate the non-linear density 
contrast. The simplest generalization does not work very well since it misses important 
features of the gravitational dynamics. These are discussed in section 4 which derives 
the main result of the paper, estimates cr^ at non-linear epochs and compares the results 
with numerical simulations. The last section discusses the implications and possibilities 
for extension. 
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2. Newtonian Gravity in expanding coordinates 

This section reviews some known aspects of clustering in a framework somewhat different 
from, though equivalent to, the standard approach. Consider a set of particles interacting 
via Newtonian gravity in an expanding universe. We confine our attention to regions with 
dimensions much smaller than the Hubble radius so that the equation of motion for the 
z*'* particle is well approximated by: 



= -E 



Gm 



(1) 



Here stands for the proper coordinate related to the comoving coordinate Xj by rj 
a(t)xj. Using this we can write the above equation as 



2d . Gm 

a ^ I X. 
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- -Xi = -V$ 
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(2) 



The right hand side of (2) has been expressed as the gradient of a scalar potential $ which 
obeys the equation 



^5(x-x0 + 



47rG Va 



a 



(3) 



Equations (2) and (3) govern the dynamics of the particles in the expanding background. 
In (3) the "mass density" is contributed by two terms, one of which is independent of x. 
This fact can be used to recast the equations in a more suggestive form by choosing a(t) 
such that it satisfies the equation : 



a AttGpo .ao _ AttG 
a ~ 3 ^~a^ " 3~ 



Pb{t). 



(4) 



This equation describes the evolution of the expansion factor in the matter dominated 
phase of a Friedman universe. Using these expressions we can write (2) and (3) in the 
form 

2d . 



Xi + 



a 



-V$; 



(5) 
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(6) 



In the above form there is explicit dependence on t due to the presence of the terms (d/o) 
and a^. It is useful to transform these equations such that the time dependence disappears. 
Introducing new dimensionless time and space coordinates r and via 



T = ln(t/T); Xj = Lqi; L^ = Gmtl/al 



(7) 
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with an arbitrary constant T, and transforming the equations, we easily find that: 

Vp = 4nJ2K^-^j)-l- (9) 

i 

This equation has several noteworthy features which we shall briefly comment upon. 

(i) The compensating negative background density, of (—2/3) in the units chosen, 
cancels the long range part of the gravity coming from the mean density. The usual 
assumption made in the analysis of Jeans instability is thus justified in these coordinates. 
One gets exponential growth in r = Int i.e., power law growth in t. More precisely, the 
analysis without the damping term would have led to modes with e'^'^ with 7^ = | . With 
this term, one gets 7^ + ^7 = f , i-^-, 1 — | or —1- The damping term has ensured that 
the growth and decay rates of the two modes are not equal and one recovers the correct 
t^/^ and t~^ factors for the growing and decaying modes. 

(ii) The damping term is also needed to recover other known results. For example, two 
bound bodies spiral inwards with a shrinking period in the q and r coordinates, this being 
just an alternate description of a bound system with a fixed size and period in r and t. It 
is also instructive to look at the free evolution of a test particle in an otherwise uniform 
universe. Writing = dqi/dr, one gets dpi/dr — — i-e., Pi oc exp(— ^r) = t~^/^. This 
is the well known result that the peculiar velocity of a particle decays as the inverse of the 
scale factor, when expressed in terms of p. 

(iii) Particles at the boundary of a void feel the strong repulsion coming from the 
locally uncompensated negative background term which pushes them further out. 

(iv) The translation invariance in r, without the damping term would just lead to a 
conventional interacting system of particles. The presence of the damping term reminds 
us that we should not be trying to apply statistical mechanics with this interaction quite 
apart from the well known difficulties associated with the singular short range behaviour 
of Newtonian gravity. 

(v) The usual calculation (multiplication by q and integrating) which leads to the 
conservation of energy for a time independent potential will now show, because of the 
damping term, that the 'energy' in these coordinates always decreases, being constant 
only in the trivial case in which all 'velocities' and 'accelerations' are zero. This decreasing 
'energy' can be viewed as the driving force behind the continued evolution with formation 
of bound structures and evacuation of voids. In fact, the form of behaviour closest to a 
steady state which is known for this system is self similar clustering (Peebles 1980). The 
various distribution functions (two particle, three particle, etc.,), describing the system 
have a particularly simple dependence on r. Writing dq/dr = p, we can have for example 
a two particle distribution function of the form 

/(gi - g2, PUP2; t) = e'^^^gie^^qi - Q2), e">i, e^^s) (10) 

Here / has been written as a function only of the coordinate difference qi—q2 after averaging 
over the mean position {qi + q2)l2. Notice that before such averaging the system is not 
self similar. 
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Finally, these equations suggest a possible way of interpreting the success of adhesion 
models (Bagla and Padmanabhan, 1993 a). This work is in progress and will be reported 
elsewhere. 

3. Connecting linear and non-linear density contrast 

We shall now describe the statistically self similar situation in more detail. Let us suppose 
that the power spectrum at some epoch has the form P{k) oc fc"', at sufficiently large 
scales which are linear. Then oc k^P{k) oc x~^'^^^^ at the early, linear stage. Since 

a grows as a in the linear theory, it follows that a\{x,a) oc a^x~^^^^\ This evolution is 
self-similar with a\ oc §"("■+3) and s = (x/t"); a = [4/3(3 + In the extreme non-linear 
case, bound structures with fixed proper radius / = a{t)x will not participate in the cosmic 
expansion. At a fixed l,cr%j^ =< {Sp/phf > must now grow as (aP/a^) = a^. The 
factor arises from the fact that the background density pb decreases as and oc p^^; 
the factor in the denominator arises from the fact that in all samples of proper radius I 
used in computing the variance, only a fraction proportional to will contribute. Thus, 
in the non-linear regime c^^(t, oc a^{t)F{a{t)x). The form of F can be determined 
by matching the linear and non- linear expressions at say, a = ac- This will lead to the 
standard result (Peebles 1980) that cr^^ oc a^[ax]~'^ with 7 = 3(n -|- 3)/(n -|- 5). Thus we 
find that 



al{a,x) oc a\t)x-^^+'^; a%Lia,x) oc a^[ax]-^ oc a^/in+5):,-Hn+3)/(n+5) (^^^ 

A graphical version of this argument is given in Fig.l. Self-similarity of this solution is 
evident from the fact that (Tnl oc s~'^ . Notice that, in the range n > —2, the non-linear 
correlation function at fixed comoving scale x grows more slowly compared to the linear 
spectrum. It follows from these relations that one can write 

= — (for cTL > ac) 

where ctc is the typical value around which the transition to non-linearity occurs. We shall 
now examine this result in greater detail. 

The above result (12) is supposed to give the non-linear density contrast at a given 
scale provided the linear density contrast at the same scale is given. Unfortunately, this 
result is of very little practical utility and suffers from several limitations as it stands. Let 
us study these limitations since they need to be tackled in any attempt to improve upon 
the relation. 

To begin with, this result is obtained by assuming a power law form for the spectrum. 
The power spectrum e.g., just after recombination, is not a pure power law in any of 
the realistic models of structure formation. Hence it is not clear what is the value of 
n which one should use. One could, of course, try to provide a local definition of n by 
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taking nefF = {din P/ dink). Unfortunately, the resulting ctjvl does not lead to even rough 
agreement with the numerical simulations. 

Secondly, notice that the actual value of depends on the value of chosen. It is 
difficult to estimate a"c from theoretical considerations with any degree of accuracy. Even 
if (7c is estimated, the above relation can be valid only for ctjvl ^ Cc and ctjvl <tC cTc- In 
the regime with aNL — cTc we do not have a reliable estimate. 

Finally, notice that the above expression attempts to relate (TNLix,a) to ctl at the 
same x and a. This is somewhat unrealistic because non-linear growth at a given scale 
cannot be expected to be a strictly local function of the linear density contrast at the same 
comoving scale. 

Let us now ask how some of the above difficulties can be circumvented by a more 
sophisticated approach. One possibility which suggests itself is to write a differential 
equation for a^L which can be generalized to the case in which the power spectrum is not 
a pure power law. This is fairly easy to do by noting that 

dlna _ 3 dlna _3{n + 3) 

dlna ~ (n+5)' ^Ina; ~ 2(n + 5) ^ ^ 

Hence, eliminating n, 

dlna dlna 2 
dlna dlnx 3 

Such an equation implies that at each epoch and scale, one is evolving a using the result of 
the self similar solution for the currently and locally applicable conditions. This equation 
is assumed to hold for a > ac- For a < ac we use the rather trivial evolution equation of 
the linear theory: 

= 1- (15) 

dlna 

These equations are to be interpreted as follows. At some very early epochs (say, at 
a = arec) we are given the density contrast a{x,arecix))- We will assume that the initial 
epoch is chosen sufficiently early so that a^ec ^ 1 at all relevant scales. With this initial 
condition the above equations can be integrated forward for all further times thereby giving 
a {a, x). This approach clearly does not require o" to be a power law. 

Because of the extreme simplicity, the evolution equation above can be integrated in 
a general form. The general solution for a > ac given by 

lnaNL{a, x) = ^ln(-^\ + F{ax) (16) 

where F is a function to be determined by matching aj^L and a^ at a = ac- Since the 
linear density contrast evolves as cr^ oc a, it follows that we can relate the ctjvl and cr^, as 
follows: 

aNL{a,x) = a^^^^{a,l); l = xa']ll{a,x) {for aNL > cTc) (17) 
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2 /3 

This shows that when cr > ctatl ^it a point x is determined by cr^ at a point I = x(jj^j^. 
In other words non-hnear evolution does introduce a degree of nonlocahty in comoving 
scale X when this approximation is used. 

While the above attempt is an improvement, it is still not very satisfactory. The 
differential equation (14) is difficult to derive directly from any quantitative dynamical 
considerations. It is merely the simplest equation that can be written down in which the 
power law index n is eliminated. That this equation is inadequate is clear from two aspects 
of our solution. Firstly, it does not make a smooth transition from linear to non-linear 
scales and secondly it does not provide us with a value for 0"^. Nevertheless the solution 
suggests a very simple ansatz for the evolution of the density contrast: It may be possible 
to take the true density contrast, crjvL(o, a;), to be a universal function of cri,[a, /(a;, ctjvl)] 

2 /3 

where /(x, (Tnl) — x for gnl -C 1 and f{x, aNL) — xcr^j^ for gnl S> 1. [In fact, a relation 
with this structure has recently been suggested by Hamilton et al. (1991) and used to 
obtain the input linear spectrum from the observed non-linear structure observed at the 
current cosmological epoch. Our emphasis in this paper is on exploring the physical origin, 
dynamical derivation and implications of such a scaling relation.] In such a more realistic 
theory, we expect / to vary smoothly from one limit to the other. It turns out that such 
a scaling relation can indeed to be obtained from the equations governing the evolution of 
density perturbations. This is task which we shall now take up. 

4. Non-linear dynamics - a better approximation 

We begin by defining the correlation function ^(x, t) as the fourier transform of the power 
spectrum: 

To an excellent approximation, ^{x^t) and (7{x,t) are related to each other by 

a\x,t)^\ [ ay,t)y^dy (19) 
Jo 

We shall now obtain an exact equation satisfied by a{a, x) by an argument which follows 
the treatment of (Peebles 1980) for the evolution of the correlation function ^. Since the 
mean number of neighbours to any particle is given by 

Nix, t) = (na^) [ 47ry^dy[l + ^(y, t)] (20) 

J o 

when n is the comoving number density, the conservation law for pairs implies 

where v denotes the mean relative velocity of pairs at scale x and epoch t. Using (19) , we 
find 
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Substituting this in (21) , we get 



2 dx 



^ ^ r 2/1 2\ 



Integrating, we find: 



(22) 



(23) 



(24) 



The integration would aUow the addition of an arbitrary function of t on the right hand 
side. We have set this function to zero so as to reproduce the correct hmiting behaviour 
(see below). It is now convenient to change the variables from t to a, thereby getting an 
equation for a^: 



d 

a— [1 + a'^{a,x)] = 



^ \^[x\l^a\a,x))] 
-ax x^ ax 



(25) 



This equation shows that the behaviour of (T^(a, x) is essentially decided by the dimension- 
less ratio h{a,x) = [v{a, x) / {—ax)] between the mean relative velocity v and the Hubble 
velocity ax — {a/a)xprop- To understand the behaviour of this equation let us consider its 
solutions in two limiting cases. In the non-linear limit, peculiar motion exactly compen- 
sates for the Hubble expansion to form bound structure; hence {v/ — ax) = 1 giving. 



(26) 



This equation has the general solution (1 + a^) = a^F{ax). In the linear limit, we know 
that — Aa'^x~^'^~^^^ ^ 1, if the power spectrum is P{k) oc k^. Ignoring compared to 
unity on the right hand side of (25), we get 



oa 



—ax 



2Aa2x-("+3) 



(27) 



Therefore, 




= (2/3)Aa2x-("+3) = {2/3)a^ oc [t-4/3(n+3)^]-(n+3) ^ ^-(n+3) (38) 



where s = {x/t") and a = [4/3{n + 3)]. The solution in the non-linear limit, (1 + cr^) = 
a'^F^ax) will be a function of s alone only if F{ax) = (ax)'* with 7 = 6/ {3a + 2); in that 
case a'^^"''x^'^ oc s-6/(3a-l-2)_ 'pj^jg recovers the 'standard results' and justifies the boundary 
conditions chosen earlier. Also note that, in the linear limit, {v/ — ax) = (2/3)cr^. This 
result in the linear regime is true for any spectrum, not for just a power law (Peebles 1980). 
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It is now clear that any general relation between gnl and ui, will translate itself into 
a similar relation for h in the two limits. Conversely, the behaviour of /i(a, x) will allow us 
to determine (7^?^, in terms of 0"^- The results obtained above suggest the hypothesis that 
the quantity /i(a, x) — [f (x, a)/ — ax] is purely a function of cr^. 

h{a,x) = h[(T'^{x)] (29) 

In the linear limit, when <^ 1, we saw that /i(a, x) = (2/3)(T^; in the extreme non-linear 
limit (cr^ ^1) bound structures would have formed in which v will balance out Hubble 
expansion; that is, v = —ax implying h = l. It seems reasonable to explore the assumption 
that h depends on a and x only through some universal function h = /i(cr^) which has the 
asymptotic behaviour derived above. Of course, this assumption has to be ultimately 
checked by comparing the results with the numerical simulations [ A word of caution on 
our notation is in order. We are using the same symbol h to denote two functional forms 
/i(a, x) and h{a'^). A more pedandic notation will be h{a, x) = H{a'^), introducing another 
symbol H. We will not do this since it should be clear from the context what we mean by 
h]. 

When h{a,x) = h[a'^{a,x)], it is possible to find a useful solution to (25). We are 
interested in the solution to this equation which reduces to the form cr^ oc a^x~^'^'^^^ for 
cr^ <^ 1. This can be obtained as follows: We rewrite the equation in the form 

- h{D)x^ = 3hD; (30) 
da ox 

where D = {1 + cr^). Suppose we can find a function F{D) such that 

F{D) = alna + pinx (31) 

with a, P, being constants. If this is true, then 

a 3 

F'dD = -da+-dx (32) 
a X 

giving a{dD/da) — [a/F') and x{dD/dx) = {P/F'). Substituting this into (30), we find 
that F{D) must satisfy the relation 



F'(r» = ^ (33) 



a-/3/^^^^^ /-(o-A) <iT^ (34) 



Or, equivalently. 



3hD J J Sh (1 + ct2) 

The solution — a^{a,x) is obtained by combining (31) and (34) We can easily find that 

n2 
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where m = (/?/«) is a constant. When /i(cr^) ~ (2/3)cr^ and cr^ -C 1, we can ignore cr^ 
compared to Ina^ and this relation gives the hnear theory result: cr|^ oc a^x^"^ allowing 
us to relate m to the index n of the power spectrum. Since we expect cr£ oc a^x"*^"^"'"^-' 
if P{k) oc /c"^, we get m = — (n + 3)/2 . This is the solution with the correct boundary 
condition which we need. It is also easy to verify that when ~ 1, cr^ » 1, we get 

oc (a^x^"^)^/^^-"^^ oc a'^iaxf'^'^^-'^^ oc a^{ax)-^ (36) 

where 7 = 3m /{m — 1) = 3(n + 3)/(n + 5). Thus the limiting forms derived earlier arise 
from the limiting behaviour of h{a'^). 

The solution we have found also exhibits the universality we noticed in the earlier 
approximations. To see this, note that the right hand side of (35) is a function of cr|(a, / = 
x{l + (T^)^/^) oc a'^x'^^ (1 + (T^)^'"/^, which is the density contrast ul{(i,1) in the linear 
theory evaluated at the scale I = x{l + a'^Y^^. Thus we can write 

-LM = e.v-j j^^^^j^^ (37) 

Since the original power law index has disappeared, this relation suggests that cr£(a,i) is 
a universal function of the correct (T^(a, x). This is precisely the kind of ansatz suggested 
by our earlier approximation. The function f{x, a) — x{l + o"^)^/^ has the limiting form, 
/ ~ xcr'^/^ for (T 1, andf ~ x for cr <g 1. 

The above relation between cr| and cr^ is indeed the correct solution to equation (30) 
even when aj^ is not a power law. This can be seen by direct substitution of the solution 
(37) into (30) or - more formally - by the foUwoing argument. Let A = Ina, X = Inx and 
D{X, A) be the solution to (30) . We define curves in the {X, A) such that 

dX 

-h{D[x,A]) (38) 



dA 



c 



i.e. the tangent to the curve at any point {X,A) has the value determined by h at that 
point. Quite clearly, along the curve, the left hand side of (30) is a total derivative allowing 
us to write 

dD ,,^,dD\ fdD dDdX\ dD , 



Integrating 



Or, Equivalently, 



I f dD 



2 [ da'^ 2 21 /.-.X 

^"^3 7 /.(.2)(i + ,2)=« <^-l\, (41) 

c 

We not only need to determine the form of the curves to fix the scale I. The equation to 
the curve can now be written as 
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dX ^ 1 dD 



giving 



3X + \nD = \n[x^{l + a^)] = constant (43) 

This shows that cr| should be evaluated at fixed I = x{l + cr^)^/^, giving us the result of 
(37). 

The actual relation between ct^vl and (Tl depends on the forms of /i(cr^). From the 
limiting behaviour discussed before, we know that /i(or^) ~ (2/3)cj^ for small cr and /i ~ 1 
for cr ^ 1. Our results for the nonlinear density contrast depends crucially on the manner 
in which h{a) reaches unity. We will expect an overdense region to expand more slowly 
compared to background, reach a maximum radius, collapse and (vivialise) to form a bound 
structure. During the collapse phase the average velocity will (in general) overshoot the 
Hubble expansion velocity. Only after varialization is complete will the pair correlation 
velocity approach the asymptotic value of (—da;). In fact numerical simulations show that 
/i(cr^) has a single maximum (and hence overshoots \v\ = dx before falling back) at about 
cr^ ~ (8 — 15) with hmax — (1-5 — 2). The simplest model for such a function with correct 
asymptotic behaviour will be: 

^/ 2^ 2 2 (1 + Aa2) 

^^^^ - r (l + (2/3)A.^) 
which has just one free parameter A. Substituting this in (37) and integrating, we find that 
alia, I) = a^{l + Aa2)(3^+2)/3(i-A)(i + ^2)-(3+2A)/3(i-A) (45) 

where — a'^{a,x). Best agreement with the numerical results is achieved for A ~ 0.36 
for which 

alia, I) = a\l + OMaY'°\^ + ^')-'-'''; / = a;(l + a^/^ (46) 

This result gives the true a'^{a,x) as a function of the linear density contrast cr|(a,Z) 
evaluated at / = x{l + cr^)^/^. 

This result does not suffer from the limitations of the earlier approximations. It gives 
the exact density contrast a^ in terms of the density contrast of the linear theory al by 
a function which makes a smooth transition from -C 1 regime to ^ 1 regime. The 
non locality of behaviour seen in our earlier approximations is also preserved in a smooth 
manner. 

The ultimate test of any such formula, of course, is based on the agreement with 
numerical simulations. It turns out that our formula does remarkably well on this count. 
To begin with, the results of the numerical simulations with scale invariant power spectra 
does show the kind of universality obtained above: The correct density contrast cr(a, x) 
is indeed expressible as a universal function of the linear density contrast aL{ci,l) where 
I = a;(l + (T^)^/^. The best fitting function to the numerical data is the multiparameter fit. 

1 + 0.0158(7^ + 0.000115(7^ ^^/^ 



^2 _ „2 

aL = a 



1 + 0.926(7^ - 0.0743f76 + 0.0156(7^ 



(47) 
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In figure 2 we have plotted this multiparameter fit as well as our result in equation (32). It 
is remarkable that our single parameter fit is reasonable for all o"^ ^ 6 and for cr^ ^ 0.2. 
The deviation in the intermediate range is due to the fact that (46) has the behaviour 
{olIoY = 1 + C'(cj^) while the actual data suggests {olIoY = 1 + 0{a^) for a\ ^ 10. 

It should be stressed that (j\ and in the above formulas are evaluated at different 
length scales. The relation of cP' to a\ at the same length scale is quite complicated and 
depends on the shape of the spectrum. For example, in the extreme nonlinear limit, 

o^{x,a) ~ 11.42/(-+5)a6/(-+5)(:,/xo)-^ ^ 11.42/(-+5)a-2("+2)/(-+5)^2 ^a;, a)(x/xo)^ (48) 

with /? = [(n + 3)(n + 2)]/(n + 5), 7 = 3(n + 3)/(n + 5) and ~ 8h-^Mpc. For a CDM 
spectrum, galactic scales have n ~ —2 and smaller scales have n ~ —3. So the nonlinear 
density contrast is a factor of about 11 larger at subgalactic scales and a factor of 11^/'^ 
higher at galactic scales. This has the effect of steepening the correlation function at the 
nonlinear end. 

5. Discussion and Conclusions 

The basic proposal of this paper is that significant progress in understanding the evolution 
of the correlation function in the non-linear phase of gravitational clustering can he made 
via a scaling hypothesis . This assumption - that the pair velocity is a function only of the 
density contrast on the same scale - leads to the explicit formulae of section 4. 

It will be interesting to test the hypothesis on a wider range of initial conditions than 
the power law cases explored so far. Given that both the velocity and the density contrast 
are averages over a variety of environments, the emergence of such a relation is not obvious 
(except in the extreme linear and non-linear regimes) and deserves dynamical validation. 
It would be interesting to see whether the form of the function h{a'^) can be obtained from 
some approximate dynamical models e.g., the Zeldovich approximation, the frozen flow 
model etc. In a widr context one would also like to investigate the conditions under which 
the universality is broken . For example, one can visualise two scenarios which start with 
the same input power spectrum, one having random phases (gaussian fluctuations) and 
the other with non-gaussian correlations. One does not, in general, expect the resulting 
correlation functions to match in the non-linear regime. Even in such situations, analysis 
in terms of deviations from 'universal' scaling might be an economical and physical way of 
compressing the information from a large ensemble of simulations. These issues are under 
investigation (Bagla and Padmanabhan, 1993b). 
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Figure Captions 



Fig. 1 A schematic diagram showing the geometrical origin of Peebles scaling rules dis- 
cussed in the text. The correlation function is plotted against the proper length in a 

logarithmic scale. The two curves indicate the function ^{ax) at two different epochs 
a = ai and at a = aie. At large scales, in which the evaluation is linear, the correlation 
function changes through two processes. The scale is stretched due to expansion and the 
linear density evolves as S oc a . At very small scales, where the evaluation is non-linear, 
the length scale does not change due to expansion, but the correlation function should 
increase as the cube of the expansion factor . Elementary geometry now shows that if the 
slope of the line at the linear end is [— (n + 3)] then the slope at the non- linear end must 



be 



3(n+3) 
(n+5) 



For a more detailed discussion see text. 



Fig. 2 The true density contrast cr^ is plotted against the linear density contrast. Note 

that the two contrasts are calculated at two different scales as discussed in the text. The 
unbroken curve is an exact fit to the numerical simulations while the broken curve gives 
the result of this paper. There is fair agreement between the analytic approximation and 
numerical simulations for density contrasts of the order of 10 and above. 
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